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THE SECOND-ORDER GRAVITATIONAL REDSHIFT 


I. INTRODUCTION 

Gravitation has always been one of the most fascinating 
phenomena of nature. Recently, a resurgence of interest in this 
field has occurred, spurred on by many new astronomical observa- 
tions. Because the force of gravity is weak compared to nuclear 
and electromagnetic forces, it is virtually impossible to per- 
form laboratory experiments involving relativistic gravitation. 
However, with advances in space technology, we have now enlarged 
our laboratory for controlled experiments to the size of the 
solar system. 

New tests of gravitation are presently being made, using 
new technology: Light bending by a gravitational field is mea- 

sured, using very long baseline interferometry to observe the 
microwave signals from quasi stellar objects, as the signals pass 
close to the sun. Powerful radars are being used to send pulses 
across the solar system to measure the time delay produced by 
the sun's gravity. Space tracking of probes flying past Mars 
and Venus have yielded valuable orbital data on the general 
relativistic corrections to the motion of objects in a gravita- 
tional field. 
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It is significant that atomic clocks have been involved in 

nearly all of these experiments. Indeed, today, atomic clocks 

offer a measuring capability within a precision of better than 
14 

1 part xn 10 , As no other instruments of comparable capability 

exist, it seems logical to further explore space-time with clocks, 
and to make measurements in terms of data expressed as time inter- 
vals. 

It appears that a direct measurement of the non-linear term 
of the gravitational field equations can now be made by the use 
of very stable clocks. A vehicle containing such a clock could 
be put into a highly eccentric orbit around the sun, or allowed 
to fall into the sun, and its frequency would be compared with 
an identical clock on earth. This is a measure of the so-called 
" second-order gravitational redshift". 

The only other presently feasible test for the second-order 
term is that which measures the perhelion advance of a planet or 
satellite; all the other current experimental tests of relativity 
can only measure first-order effects. This frequency shift 
experiment would be a distinct,^ additional test for the nonlinear 
term. A great advantage of this experiment is that it does not 
depend at all on the integrated characteristics of the orbit. 

As we will show, the only requirements are the knowledge of the 
instantaneous gravitational potential and the velocities of the 
clocks. 
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At present, a highly accurate test of the first-order gravi- 
tational redshift is in development, using a terrestrial rocket 
probe. Some of the techniques perfected for that experiment will 
be useful in a heliocentric second-order redshift experiment. 

The terrestrial probe experiment uses atomic hydrogen maser 
clocks in the probe and on the ground that are connected by 
microwave signals. The nongravitational effects due to the rela- 
tive motion of the probe and ground clocks are accounted for by 
continuously measuring the path length in terms of the phase of 
coherently returned signals originally transmitted to the probe. 
The phase variations due to path length changes are then auto- 
matically removed from the clock signal. 

The terrest^rial experiment is expected to test the so-called 
principle of equiv'^alence. This principle, first stated by 
Einstein in 1907, asserts that all freely falling, nonrotating, 

t 

"sufficiently small" laboratories are equivalent; i.e., all the 
laws of physics will appear the same in all freely falling, infin- 
itesimal laboratories. At present, the best tests of the princi- 
ple of equivalence are the elegant experiments performed by 

R. V. Pound and co-workers, using Mbssbauer emission and absorp- 
. 57 

tion from Fe . These experiments, performed over a vertical 
distance of 75 feet, have verified Einstein's postulate of equi- 
valence to 1 %. The terrestrial probe experiment seeks to extend 
the distance to 10,000 miles, and the accuracy to 20 ppm. The 
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successful completion of this test will do much to put the equi- 
valence principle on a sound experimental basis, and establish 
the validity of the cornerstone of Einstein's General Theory of 
Relativity. 

However, the terrestrial probe experiment will not be suffi- 
^ ciently sensitive to permit distinguishing between the several 
theories of gravitation now in contention. Other presently viable 
theories of gravitation have the principle of equivalence as a 
basis. The first-order gravitational redshift can be derived 
immediately from the principle of equivalence, alone; a knowledge 
of the field equations of the theory is not necessary. Thus, it 
is. important to emphasize that a measurement of the second-order 
redshift, the nonlinear term, is a test' of the actual field equa- 
tions of a theory of gravitation, and could help establish the 
validity of such a theory. 
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II. FREQUENCY SHIFT IN A GRAVITATIONAL FIELD. 


We derive the exact expression for the gravitational frequency 
shift of light emitted by a source at one position .and received 
by an observer elsewhere. The rigorous, covariant approach is 
used in order to avoid any ambiguities in the definition or inter- 
pretation of any terms appearing in the final expression. 

The ratio of the frequency of- light emitted by a source (s) 
and received by an observer (o) can be given as: 


V, 


V ' 
. 0 


(P • v) 
(p • v) 


'Pg ''“>0 


( 1 ) 


dx 


a 


a 


II ct _ dx , 

where p^ = xs the photon "four-momentum, v = xs the 

"four-velocity" of the source or observer, and ds is the invariant 

1 

xnfxnxtesxmal element of length xn four-space. ds = cdx , where 
is the so-called "proper time" for a particle; for light, 
ds = 0, and, in that case, ds is simply an arbitrary variable. 

The symbol ot represents numbers from 1 t o 4. Similar subscripts 
and superscripts imply summation over that index; this is the 
usual convention. 

The proof of Eq. (1) is quite direct. (P . • is the product 

of two four-vectors; it is a scalar — an invariant quantity. Being 
invariant, it has the same value in any coordinate system. We can 
easily find the invariant value of (p * v) by considering an 
inertial system in which the particle is at rest. (By the principle 
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of equivalence, one can always make a transformation into an iner- 
tial system at any point in a gravitational field.) In this sys- 
tem, from special relativity, we have, with c = 1, P,v • " 

^‘Lj.ght 

= (0, 0, 0, 1), so that (p » v) 


(Pi ' P^> E = hv ) and v 


a 


Particle 


(P^j v“) - hv. Thus, 


(p • v) 

(p • V) 


hv 

c 

ET 


V 


Eq. (1) is quite general; we now apply it to the special case 
of the spherically symmetric gravitational field. The general form 
for this field, in the (r, 6, ‘P, t) coordinate system, is given as^ 


ds^ = B{r)dt^ - A(r)dr^ - r^{d<j)^ + s.in^0d0^) 


( 2 ) 


Upon solving the geodesic equations of motion, for either light or 
a particle, it is found that the motion can be confined to a 
plane, (p = constant (chosen as = it /2), just as in the Newtonian 
case of central force motion. Eq. (2) can then be taken as 


ds^ = B(r)dt^ - A(r)dr^ - r^d0^ 


(3) 


Evaluation of p_ . , . and v . . 
^Light PartiOle 


The geodesic equations for 6 and t give 


ds 


B 

QS 


= q = constant 


= k = constant 


(4) 

(5) 
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If, for a moment, we consider a test particle at infinity, of 
finite mass m, with velocity v as measured af infinity, Eqs. (4) 


and (5) can be written as 


2 de 

r — = = 

as ^ 


£v 


/ 1 - V' 


B II + ^ = k = 

ds as 


( 6 ) 

(7) 


/ 1 - 


V 


where ^ is the classical impact parameter of the particle as 
measured at infinity. Therefore, for light, where v = 1, we have 


2L = ii 
k ^ 


(8) 


Using Eqs. (3), (4), (5), (8) and ds = o for light, we have for 


dx 


a 


^Light ■ ds 


p" = H 


p' “ ^ 


k 

B 




2 = ^ 

P ds 2 


/i - 


dx 


a 


and, finally, for ^ 


Po = goo P° =- ^ 

Pi = Ull P^ = - k y ^ 

P^ = 92 2 P^ = - . 


/F /TT 




(9a) 

(9b) 

(9c) 


For a particle, ds o. Using Eq. (3) , we find for 

V = dx . 

Particle ds 


a 
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,0 _ 


dt 

ds 


/ 


B - A 


( 10 a) 


1 _ 


V 


2 _ 



dr 

dt 

ds 

dt 

ds 


dQ 

dt 

ds 

dt 

ds 


V V 
r 


Vq V 


(10b) 

( 10 c) 


dr 

dt 


d 0 

where v e ^ and. v g = . Therefore, 


-y a 

p • V = p V 

a 


= k (1, - 


T~ f- 

B / 1 - 


t-2 / 


- £) v° (1,V^.,V ) 


= k 


/V- 


-A 2 2.2 

B-Av^ -r V 0 . 


-'T A” 


2il v - l.'ii) 

2 r 6 


In the absence of gravity, A = B = 1, and the ratio of the frequency 


sent by a source, moving with velocity (v^, V 0 > at a point r, versus 
the frequency received by an observer at rest is 


l-/l-— 


(p*v) 
(p *v) 


2 Vr~- (rvg)] 




2 2 2 
'^r ^6 


1 

k 


[1 


- /TV 


- - (rvg) 


/ 1 - 


7 7 2 

V ^ - r'^Vn 

r 6 


This is just the special relativistic doppler equation, as we 

should expect: The numerator is (1 = v^ , cos 0), as is 

^ Particle 

verified, by a simple geometrical construction. The denominator 
is /l - v^ , the inverse Lorentz factor (c = 1) . 
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Eq. (11) shows that the presence of a gravitational field 
"modifies" the well-known doppler-shift equation through the func- 
tions A(r) and B(r). The parameter £ is the single constant of 
the motion for light. In the absence of gravity, it is clearly 
the impact parameter of the photon; in the presence of a gravita- 
tional field, this unambiguous physical interpretation of £ can 
only be made at infinity, as will be discussed below (Section VII) . 
No physical interpretation can be made for r “, even though £ 
is still a valid mathematical constant of the motion for the 
photon. 

We can write Eq. (11) in a more compact form, if we define 


T 


/ 


B - A V ^ - r^v„^- 

i O 


t = 




( 12 ) 

(13) 


(c will always be taken = 1, unless otherwise specified), and 

C = ( /T/ 1 - , i , 0) . 


(14) 


Eq. (11) becomes 


•V = kyd"?'?) 


(15) 


In the absence of gravitation, y reduces to the Lorentz 
factor, and e is the unit direction vector of the photon at the 
point r. We emphasize that, when A and B ^ 1, y and e are 
modified by these functions, and such a simple interpretation 
cannot be made. 
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III. THE DOPPLER-CANCELLING TECPMIQUE 


Because of the motion of the space vehicle relative to earth, 

the first-order (in v/c) doppler effect is clearly -the largest 

contribution to the signal sent by the vehicle clock. The lowest- 

2 2 ~ GM 

order gravitational contribution is of order v/c ( ~ ). 

c^r 

It is thus essential to set up an experiment that eliminates the 

10 , 11,12 

first-order terms. This can be done by placing both a 

clock and a transponder on the probe. The transponder serves as 
a reflector that allows us to measure the two-way distance in terms 
of signal phase. By combining the received signals appropriately, 
we can remove the first-order doppler terms, as follows: 

We first consider the situation where the earth station is 
equipped with a clock and the probe is equipped with both a clock 
and a transponder. The earth station sends out a signal of proper 
frequency v at proper time t^ (with the station at position 1, 
with velocity ) to the probe's transponder, which receives it 
at time t 2 ( position 2, velocity ^2 ). The probe transponder 
re-transmits the signal to the earth, along with the probe's clock 
signal; both signals are received at the earth station at time t^ 
(pos-ition 3, velocity Ba ). (Figure 1) 

Eqs. (1) and (15) give the frequency shift of the clock signal 
received from the probe as 


^ 3 
V2 


Clock 


Y 3 ( 1 -^ 3 * £ 2 a ) 

Y 2 ( 1 -^ 2 ’£23) 


(16) 
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where is the "direction vector" of the light signal from 

position i to j. (p 2 • £23 is evaluated at position s, P3 * £23 
is evaluated at position 3, etc.)* 

In a similar manner, the frequency shift of the earth-probe- 
earth transponder signal, received on the earth at time t^ is 


V3 

VT 


Transponder 


V2 V3 


Y3 

Y7 


(1_^2*£i 2 ) 
(1-^1 *ei 2) 


( 1-^3 *£ 23 ) 


(1 


2 3) 


(17) 


The doppler-cancelling technique consists of combining these 


two received signals in the form: 
Ay 

V 



1 V 3 

= 

Earth V 2 

Clock 2 Vi 


Transponder 


1 

2 


(18) 


The observed — is very small, and contains the information we 
seek; it contains only terms in 8 and higher order. 
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IV. THE DUAL TRANSPONDER SYSTEM 


It is also possible to equip the earth station with a trans- 
ponder, and use the doppler-cancelling scheme for earth-based 
clock and transponder signals received by the probe . (This infor- 

i 

mation is subsequently transmitted to the earth station for analy- 
sis.) This dual transponder-clock arrangement provides twice the 
number of observables for a self-consistent analysis of the experi- 
ment; the added information gained from the second system permits 
a more direct and tractable means for explicitly evaluating the 

various parameters in the experiment. The dual system would also 

12 

aid in eliminating possible nonreciprocal propagation effects. 

For the probe-based system, with a signal sent from the probe 
at. time t^ (position 0, velocity ), received at earth at time 
t^, and received back at the probe at time t^ along' with the 
earth-clock signal, the doppler-cancelling technique gives 
Av 



_ 1 Vz 

* 

Probe Vi 

Clock 2 Vo 


Transponder 


1 

2 


(19) 
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V. THE SECOND-ORDER REDSHIFT IN A SPHERICALLY 
SYMMETRIC GRAVITATIONAL FIELD 


We will consider the gravitational field of the sun as spher- 
ically syimnetric in the first approximation. A possibly signifi- 
cant correction to this model will be considered later. The solu- 
tion of the field equations of general relativity for the spheri- 

2 

cally symmetric case was first given by Schwarzschild shortly 
after Einstein's publication of the basic theory. In a plane ( <)> = 
the exact solution is given as 

ds = B(r)dt^ - A(r)dr^ - r^d6^ , (3) 

with 


B(r) 

= 1 _ 2GM 

r 

(20) 

A(r) 

= 1/B(r) 

(21) 


3 4 

Eddington and Robertson put the spherically symmetric solu- 
tion in a more general, parameterized form in order to show more 
explicitly the contribution of the various-order terms in (GM/r) . 
In the "standard" (r, 9, t) coordinate system, we have 

ds^ = B dt^ - A dr^ - r^d0^ , 


/ 2 ), 
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with 


A = 1 - 2a — + 2(g-aY) + 

^ V.2 


B = 1 + 2y + . . . 


( 22 ) 

(23) 


where ot, 3^ and Y are unknown, dimensionlebs parameters. 

If the Einstein field equations give the physically correct 
description of the spherically symmetric gravitational field, then 
the solution would be the exact Schwarzschild .solution, and in the 
parameterized representation we would have a = 3 = Y =1. This, 

of course, is what is to be tested. In contrast, the Brans-Dicke 

5 D wfi 

theory would have a = 3 =1, and Y = “ , where w is the 

o)+2 

unknown dimensionless parameter of that theory, which governs the 
admixture of the scalar field to the tensor field. 

Actually, a = 1 simply follows from the empirical definition 
of the mass, M. We will thus set ^ = 1 in all further calcula- 

tions. 

Earth-Based Station 

Combining Eqs. (16) - (18), we have for the doppler-cancelled 
signal at the earth station: 

Hj. 

Transponder 


V 


= ^ 

1 ^ 

Earth Vz 

Clock 2 Vi 


1 

2 


= li. (i~Ae ' ) _ ^ Xl ( 1-Ap ' ) (1-Ae * ) _ 1 
Y 2 (1-Ap) 2 Yi (l“Ae) (1-Ap) 2 

_ (1-Ae* ) r Y3 _ 1 Yj. ( 1-Ap ' ) 1 _ 1 
(1-Ap) '^Y 2 2 Yi (l-Ae) ^ 2 , 


C24) 
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, and 


where Ae ' = P3*£23 , Ap = p 2*£23 , Ae = pi*£i 2 

Ap'=^2*£l2 

We will initially assume that the directions and vector velo- 
cities do not change during the signal propagation time (i,e., 
infinite propagation velocity). In this cas^, 33=3i,Y3 =Yi , 
£2 3 = -El 2 , which gives Ae' = -Ae and Ap' = -Ap, Eq'. (24) becomes 


Av 

V 


Earth 


(1+Ae) r Yi 1 d-Ap') T 1 

d+Ap’ ) *• Y 2 2 (1-Ae) ^ " 2 


(25) 


Expanding to fourth order in 3 (or, equivalently, second 
order in Gl^r) , and using the Eddington-Robertson parameterized 
form for the gravitational field, we have (Appendix A) : 


Av 

V 


Earth 3 2 )]+ Ae(Ap'-Ae) 


- (Ap'-Ae) [( 


m _ m_ 
ri “ rz 


(26a) 

)• + I (3i^- 32^) + Ae(Ap'-Ae)] (26b) 


+ m [ Y + (iii+Bii)] ^ 26 c) 

‘ ' ri rz 2 ri 2 r2 2 ri rz 


+ i [ 3 3 i“-32‘‘-23i^32^ ] 


(26d) 


+ [ (6-y) (- 


3^v ^ 3 _ 1 _ _1 

2 2 22 2’2 ^ 1^2 
r 2 ri ri r 2 


] (26e) 


+ (Ap'-Ae) { Ap' [(^- + T (Bi^-32^) 

• r 1 r 2 2 


+ Ae (Ap'-Ae)] + Ae^ } 


(26f) 
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where G has been set = 1 (GM/c^ m) and ^ ^6 ^ 

i i 

Eq. (26a) is the second-order term, (26b) is the third-order term, 
and the rest are all fourth-order in ^ (or second order in m/r) . 
The nonlinear Eddington-Robertson parameter, g , characterizing 
the second-order (in m/r) purely gravitational redshift, appears 
only , in the term (26e) . 

Probe-Based Station 

A similar analysis can be made for the probe-based doppler- 
cancelling system. Eq. (19) becomes 


Av _ d-Ap') , Y 2 1 Y 2 d-Ae"), 1 

Probe (1-Ae) ^ Yi 2 Yo (1-Ap")-' 2 " 


(27) 


where Ae" = and Ap" = ^o*£oi . For an infinite 

signal-propagation velocity, ^2 = l^o , Y 2 =Yo ,eoj=-ei 2 , 
which gives Ae" = -Ae and Ap" = -Ap' . Eq. (27) then becomes 


^ = (1-Ap* ) r li _ 1 (1+Ae) , _ 1 

V Probe ( 1-Ae ) ^ Yi 2 ( 1+Ap ' j 2 


(28) 


Eq. (28) has the same form as the earth station relation, 

Eq. (25), with the substitutions Ae (-Ap'), Ap' (-Ae) and 

the (1, 2) indices reversed. The fourth-order expansion of 

Eq. (28) can then be found from Eq. (26) using these substitutions 
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iproti© 


. ^ ^ I^p 

X v-y (29^> 

^ ■?L-^ + 5 

1 


+ 62^3, — "^ri 


r 2 


( 29 ^) 


1 t 33i X X ^ ^ 

+ '8 ^ ^ 


•9 2 

•9 2 ^ 

4> -v" rt 


.■\ (39 e) 


_j^ ' ^ 2' ^ 

^ TCi ^ 

, < ?,e U ?7 “ 



n 



VI. FINITE SIGNAL TRANSIT TIME 

We will, in fact, have a finite go-return time for the propa- 
gation of the signal. For a heliocentric probe, there would only 
be a relatively small displacement of the earth or probe during 

t 

this time, as compared with the total displacement over a whole 

-5 

period of its motion (approximately 3 parts in 10 for any of the 
vector positions and velocities, at best). Therefore, r, ^ , and 
Efj would not change appreciably in direction during the signal 
transit time. We will use a perturbation analysis to investigate 
the additional contributions due to the finite displacement of 
the earth or probe during the signal transit time. 

We will assume here that the earth moves in a circular orbit 
around the sun and ignore the earth-moon motion and earth rotation. 
We write' 

ti = ti + t , (30) 

and = l3i|^ this case. Defining e by 

e23=-ei2+E , (31) 

we have 
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Ae ' = 1 3 • e 2 3 

= (^1+^) • (-£i2+e) 

= - ^i*ei2 + - £i 2*^ + e*^ 

= - Ae + ^i*e - £i 2 *^ + £•? / (32) 

4 

and 

Ap E B^ 2*£2 3 

= ^2 • (-e 1 2+e) 

/ 

= - ^2 •£ 1 2 + ^2 *e 

= - Ap' + ^ 2 *e . (33) 

All the terms have been retained in. these expansions, since we are 
carrying our results through to fourth order in 6 . 

We substitute Eqs. (32) and (33) into the general relation 
for the earth station, Eq. (24) . Further calculations (Appendix B) 
ultimately yield 

^ Earth = ' =■ '='5- <26)J (34a) 

+ I [(£ 12 *^) + (^ 2 -^i)'e] (34b) 

- j [(si 2 *^)Ae + (B 2 *e) Ap' - (^i*e)Ae + (e*^)] (34c) 

■+ [(ei 2 -^) + (^ 2 -^i) -e] [ (^ - ^) + y ( 61 ^- 32 ")] ~(34d) 

r 1 r 2 / 

+ ^ ^ ~ (3i*e) [Ae(2 Ap'-Ae) + 
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(34e) 


+ (32*e)[2 Ae (Ap'-Ae) + Ap'^] 

+ (ei 2 *'^) [Ae (2 Ap'-Ae) + (e-t) Ae 
+ (^2*e) (£i 2*6) + (32*e) (32-3i)*e} . 


Eq. (34b) is of lowest order 3^ ; (34c) is of lowest order 3^ ; 

(34d) and (34e) are of order 3** . This is relatively simple to 
verify: 


To find the order of 6 , we can rewrite t, = t: ->■ t as 

P 3 ~ A T, where A T is the signal go-return time. For 

the earth, moving with velocity v in a circular orbit around the 

■ 0 ^ 1 2 

sun, with "radius" R, ~ . With AT “ q » we have 

-V 2 -2 ' ■ ' 

t ~ • A T a i ^ ^ ^ * This simple analysis 

c 

shows that the lowest order dependence for <5 is C?^(3^ ). 

To see the dependence ofe =£12 +£ 23 , we will use a 

6 

plane geometry argument ; this is not rigorously correct (as dis- 
cussed in Section VII) , but it will at least give the lowest-order 
dependence. We have 


e = 


£ 1 2 


£23 
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with r 


r + c A T ? 

1 1 


e 


(r 2 -r i) 


~ (r 2 -r i) 


- cAt^i 


cAT3i •£! 2 
! r-2-^r 1 


Expanding the denominator of the second term yields 


~ cAtB^i 
e ~ 


r 2 -r 1 


- e 1 2 {■ 


:At^i * £ i 


-) 


Ir 2 -r 1 

With A T ~ 2 I I /c* we finally have 

e ~ 2^1 - 2 £ i 2 (^ i *£ i 2) 


e is at least of order B . 

We see that the inclusion of the finite signal-propagation 
time introduces further terms of all orders in the doppler- 
cancelling data; these terms must, of course, be carefully treated 
on the same footing as the pure doppler and relativistic effects 
calculated in the preceding section. 

A similar analysis has been made for the probe-based station. 
In this case, no restriction to circular motion can be assumed 
for the probe. With r 2 =:?o+P/ tz = to + , ei 2 =-eoi 

+ e' , Appendix C gives 


V 


= (— ) 

Probe V o 


[ = Eg. (29)] 

- i [(ei2-^') + (l2-3i)-E'] 


(35a) 

(35b) 


+ i { [3 Ap'-2 Ae][(e, 2 *^') + (^ 2 -e')l 


r • p 

+ (Ae-2 Ap')(^i*e') - Te"<S') + ^ 

r 

0 


- (^u •?' ) } 


(35c) 
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2" ^ •?' 




*P' + (Jj-?M ,5 

^ + (?2-? , - P'<^P'-Ae, 


Rl 

Jr~* f / — P . o ^ 

;TT - I (£i:pj 2 1 

flo4« 


^ Ue-^Ap 


'> + 2 


Ap'^j 


2 (Ap ' ~^Q ^ ^ 0 * Q 


J } 


^0^ 


)] 


(35d) 


f3Sb) is Of , 

°f iowsst-ofdsr fij 

‘*- i3 „f „ . '■ <35 c) i3 ^ 

r * * S', (p , ^ ^“''^“t-ordsr ga'. 

^» + OA J. j ^ ^ Of re ^ 

’^0 + 2 \ i * , ^ ' 'or rj = 

' ^“ - ra I, j 
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VII. THE REALITY AND INTERPRETATION OF COORDINATES 
IN THE SECCND-ORDER REDSHIFT EXPERIMENT 

Throughout this work, we have used the so-called "standard 
(r, 6, <j), t) coordinate system". While this or any other coordi- 

nate system is quite valid to work in, within the framework of 
general relativity, there is the question of whether the (r, 0, <j), t) 
system physically corresponds with our familiar Euclidean concept 
of a spherical coordinate system. In fact, it does not, and the 
lack of correspondence with our familiar geometrical concepts must 
be accounted for carefully in the analysis of the second-order 
redshift experiment. These questions must be carefully considered 
when we analyze the tracking of the probe from the earth. Posi- 
tions, velocities, angles, etc., will all be affected. 

We can illustrate the nature of the problem. Let us consider, 
for a moment, the exact Schwarzschild solution in a plane, 

ds^ = (1 - dt^ - (1 - dr^ - r^de^ . 

The infinitesimal elements of spatial distance in the r and 
direction are given by^ 

1 

= (1 - ^ dr (36) 

dS,g = r de . (37) 

These expressions are found by setting up an operational defini- 
tion for the measurement of spatial displacement at a point in a 
gravitational field, using light signals.^ For the static 
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Schwarzschild solution, these d it happen to be simply the spatial 
components of Eq. (3), but this is not generally true for an 
arbitrary gravitational field. 

Using Eqs. (36) and (37), we can find the ratio of the cir- 
cumference of a "circle" at coordinate r from the origin, to the 

radial distance. This is 

2 r 2 tt r 

Jditg / /dit^ = r.j^de / ^ 

0 0 0 r 0 (1-— )2 ' 

2T 

= 2rr T 

2rr\ — 

0 ( 1-— )2 
r 

The integral in the denominator is obviously not going to be equal 

to r, but some value greater than r (since 2m < r) . Therefore, 

the ratio of the circumference to the radial distance will be less 

than 2 7 T / in disagreement with classical Euclidean geometry. The 

"radial coordinate, r" cannot be interpreted simply as the straight 

line distance from the origin to a point r. 

7 

Further study shows that there is indeed no way physically 
to interpret the coordinates in terms of our familiar geometrical 
concepts. The coordinates are a valid means of cataloging the 
points and events in space-time, but they can only be given a 
real physical meaning if they are somehow eventually written 
completely in terms of observables (proper- clock time intervals 
and proper frequencies, in this case) . This is the self-consistent 
analysis of the redshift experiment that ultimately must be made. 
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We now. consider the specific equations of the second-order 
redshift experiment. In the compact notation defined by Eqs. (12) - 
(14), the frequency was written as 


E p' 


= k 


(1 -t‘ 


However, it is incorrect to say that e is simply the plane direc- 
tion cosine of the photon, since its explicit foirm in the (r, 0, <l>, t) 
"standard" coordinate system is 



Only when A = B = 1 will e be the actual classical direction cosine 
for the photon. In the presence of a gravitational field, where 
A and B ^ 1, the form of e is modified, and no unambiguous geo- 
metric association can be made. 


Indeed, even £ cannot be associated with the classical impact 
parameter of the photon, when a gravitational field is present. 

By combining Eqs. (3), (8), (20), and (21), we find the r - 0 

equation for light (ds ^ = o) in the Schwarzschild field to be 



(— )" 



The impact parameter r. 


is the point where 


d 0 


we have 



or 



(38) 

0. From Eq. (38) , 


( 39 ) 



7 ^ r except as r “ . We see that ^ cannot be considered the 

"impact parameter", except at infinity; it is merely a constant 
of the motion for the photon, and no deeper meaning can be attri- 
buted to it. 

All the parameters appearing in the redshift equations 
( £, r., g., etc.) must eventually be evaluated rigorously in 

terms of the observables, using the redshift and light propagation 
equations in a self-consistent scheme. 

We note that in the above discussions of the various-order 
terms appearing in the doppler-cancelled data, we have assumed 
the . terms to be only of order 3.. There are, of 

course, higher order corrections to , involving terms 

in m/r, due to the presence of A (r) and B (r) ; we have not 
explicitly expanded them out in this present work, in order to 
avoid further complicating the equations. These questions will 
be addressed in the tracking study phase of this work. 
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VIII. FURTHER CORRECTIONS TO THE REDSHIFT ; 
SOLAR QUADRUPOLE MOMENT 


There are many corrections that should be included in any 

realistic analysis of the second-order redshift experiment. Of 

particular importance is the possibility of a solar quadrupole 

8 

moment, as suggested by Dicke and Goldenberg . We will now con- 
sider this correction, reserving such effects as the earth's rota- 
tion, etc. , for a later analysis. We will use the "weak-field" 
approach to find the quadrupole contribution.. 

In a"weak-f ield" , the general metric components for the gravi- 
tational field can be written as the classical, special relativistic 
term (- the Kronecker delta function) , plus a small perturba- 

9 

tion, Y / i*e., g = -6 + y .It was shown by Einstein , 

’yv yv yv 'yv ^ ' 

that, if we define 


Y 


yv 


Y ' 

' y V 


-la 


yv 


Y' 


PP 


(40) 


then 


I / 4- \ 4G 

Yyv' (x,y,z,t) = - — 

o 


Tyv (x* ,y’ ,z' ,t-R/c) 
R 


dv' 


(41) 


where is the energy-momentum tensor, with the primes denoting 


the source coordinates, and R = [ (x - x ') “^ + (y - y ') ^ + (z - z*^)^]^. 
To lowest order in v/c, we have 


yv 


~ T 


yv _ 


= , P 


dx ^ dx 


V 


ds ds 


(42) 
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where p (r') is the density. For a nonrotating mass distribution, 
we have 


T = - p (r ' ) 

H 4 


(43) 


as the only nonzero component. In this case, Eq. (52) shows that 

the only nonzero component of y ' will be v 

pv 4 4 

We combine Eqs. (41) and (42), and expand ^/R in spherical 
coordinates: 


I _ 


4G 


4 4 


' ■ , 2/ 

P(r')[ 7 ^ (^) Pn(cos x) ] dv' , 

^ ii=o 


(44) 


where x is the angle between a source point (r', 6', ({)' ) and the 

observer point (r, 6, cf)) , with 


cos X cos e cos 9* + sin 0 sin 0’ cos ((}>-(})’) . 


(45) 


Expanding Eq. (44) , we have 


Y 


4 4 


4G 

c^r 


y" p(r') [ 1 + (i) r* 


cos Y 


+ (^j Y *^'^(3 cos^y-1) + ...] r’^sin0* dr'd0'd4>' 


Inserting Eq. (45)‘ into (46) and integrating over 0 (o '.it ) , 
(o — *■ 2it ), and r, we find 


Y 


4 4 


4GM 

2 

c r 


2G 


[ I (1 - ^) + I (1 - ^) + I (1 
c-2r3 X ^2 y ^2 z 

1 2P 

+ . [ xy P + xz P „ + yz P ] , 

2 5 xy xz ^ yz ^ 

c r 


3z ' 


•) ] 


( 47 ) 


. (46 
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where m= P(r') dv' , I. = P(r') x.' dv' , and P.. = p<r')x.'x 

1 1 1 ' 

If we consider the sun to be an oblate spheroid around the 

z axis, then I = I and P. . = 0. The nonzero components of q. . 

X ■ Y ID ^iD 

are then, to lowest order. 


1 1 


’ 4 4 


= 1 - 


2 2 ^33 

2GM 


= - 1 - 


2GM 

2 

c r 


2 ^ h <1 - • 

c r c r r 


(48) 

(49) 


Confining our attention to the plane z = 0 ( <J) = tt /2) , g . becomes 

4 4 


_ , 2Gfl G ■ ^ T \ 1 

^44 ^ “ 2 2 ^Z 3 • 

c r c r 


(50) 


8 _ -5 3 

According to Dicke and Goldenberg , I - I ~ 4 x 10 M©R© , 

X z 

and Eq. (50) can be written as 


4 4 


~ 1 - ^ - (4 X 10"^) ~ (p®) 


(51) 


where our previous notation (G = C = 1) has now been used. 

We have carried out this analysis only to first-order in m/r, 

as this is all that is necessary; even if Dicke and Goldenberg are 

correct, the quadrupole term in Eq. (51) would only be of the 
-15 -1^ 

order of ~10 - 10 ; i.e., of order 3"* . The possibility of 

such a quadrupole contribution should be included in any analysis 
of the measured data. 


. ' dv ' . 
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IX. MODEL ORBIT CALCULATIONS 


We have used a computer to calculate the time dependence and 
magnitude of all the terms in the second-order redshift experiment 
(with Y =3 = 1) , for various orbits of a heliocentric probe, 

such as the proposed NASA-ESRO heliocentric satellite mission 
(Figures 2 - 4) . These calculations have been made for an earth- 
based station that moves in a circular orbit of radius 1 AU around 
the sun. Even though this model is idealized, as are the assumed 
probe orbits, we can still obtain a reasonable estimate of the 
magnitude of the various terms as the mission progresses. 

These calculations show that the combination of sufficiently 
accurate and stable clocks, and judiciously chosen orbits, make 
a second-order redshift experiment quite feasible. 
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APPENDIX A 


Expanding 


V 


i Earth 


( 1 +Ae) 

d+Ap') 


Yi 

yT 


1 

2 


( 1 -Ap* ) 
(l-Ae) 


1-1 


(Al) 


(1 + Ap 


and (1 - Ae) to fourth-order in 3 , we have 

= ( 1 +Ae) (l-Ap'+Ap' ^-Ap' ^+Ap' **) X 

Y 1 1 1 

[ — - i d-Ap') d+Ae+Ae^+Ae^+Ae'*) ] - 4 
Y 2 ^ ^ 


= ( 1 -Ap ' +Ap ' ^”Ap * ^+Ap ' **+Ae-AeAp ' +AeAp * ^ 


-AeAp ' ^ ) X 

[(Yi/Y 2 “ 1 ) - (-l+Ae+Ae^+Ae^+Ae** 

-Ap' -Ap ' Ae-Ap' Ae ^-Ap ' Ae ® ) ] - ■j 


= ( 1 -Ap ' +Ap ' ^-Ap' ^+Ap ' ‘*+Ae-AeAp ' +AeAp ' ^-AeAp ' ^ ) x 
(Yi/t2-l) 

+Ae (Ap ' -Ae) -Ae (Ap ' -Ae) ^+Ap ' Ae (Ap ' -Ae) ^ 
+AeMAp'-Ae) . (A 2 ) 
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Explicitly, from Eqs. (12), (22), and (23), is given as 


2m 


Yl = [1 - ^ + (6-Y) ^ - Br " - f 

i i i i 


^i 


1 

2 


n 2m , /o . 2m^ „ 2 

[1 - — + (B“Y) — -.3. •* 

""i r_. 2 ^ 


1 

2 " y 
r . r , 1 ^ 

1 1 ^ 


2m „ 

Y — V 


(A3) 


where 3 • ^ = 3 ^ + 3n (The unsubscripted 3 and y the 

^i ^i 

Eddington-Robertson coefficients . ). 

Expanding to second order in m/r yields 


Yi 

yT 


- 1 


1 ,2m \ 2m^ . « 2 . 2m „ 

I (8-y) + 3i^ + Y :r- 3 


2 'r 


1 


ri 

2m ^ 


|H + (0-y) 


r, -r, 


I ^ _ Y ft 2 1 

2 ra ra 


. 1 fl2m^ , 12m ft 2 , IQ 4 

+ 8 t-l~ ^ FT ^ ■ 

ri rz 

2m 


4 4m^ _ in ft 2 

7 ~ rz 


T [” + ^ ^ 3 - 2 " + 3i^32^J 

^ ra rj 


32**) 


3^/ 3 '^^ 

+ m [y(-_-™) 


1 ^ 
2 ra 


T 32^ Br 

2 ^FT ^ 


+ ^ ~ 3a** “ 


+ [ (3“y) ( 


^ 3 JL_ ^ 1 1 

y . z ^,2 2 a 2 a rira 

3: 1 r 2 . 


3. (A4) 


32 



The first line of Eq.- (A4) is of second-order in 3 ; the other 
lines are all fourth-order, with the second line being a mixed 
doppler-gravitational term, the third line a pure doppler term, 
and the fourth line being the purely second-order (in m/r) redshift 
term. 

Combining Eqs. (A4) and (A2) , we finally obtain 


Av 

V 


Earth 


(?- - T ( 3 i^- 32 ^) + Ae (Ap'-Ae) 

XI X 2 ^ 


+ 

- (Ap'-Ae) [(^- 7"^) + i ( 3 i^- 32 ^) + Ae (Ap'-Ae)] 
r 1 r 2 2 


^ 3, 


+ m [y(- 


r 1 


3 


3 


.3, ^ . 3 ? 


_ ^ ^ _ ±. ^ .. i ( . 2 4. 1 \ 1 

ri rz ' 2 ,ri 2 Vz . 2 'r ^ ^2 ^ 


+ ^ [33i'*-32‘*-23i^32^] 


l\i3 1 1 1 1 

+ m [ (3 -y) ( ) + 9 - - - 5 - ■ 

2 ^ _ 2 ^ . 2 ^ ^ 2 3: 1 r 2 


r2 r 1 


ri 


r 2 


+ (Ap'-Ae) { Ap' [. (^ " + 7 ( 3 i^- 32 ^) + Ae (Ap'-Ae)] 


+ Ae^> . 
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APPENDIX B 


The general relation for the doppler-cancelled signals received 


at the earth station is 


Av 

~ 


Earth 


(1-Ae* ) r ^ i'll d-Ap') 

(1-Ap) Y 2 2 Yi (1-Ae) 


From Eqs. (32) and (33) , 


1 

2 


Ae' = - Ae + 3i*e “ 


r 


Ap = - Ap* + ^ •£ 

2 


(Bl) 


(B2) 

(B3) 


Substituting Eqs. (B2) and (B3) into Eq. (Bl) , we have 


Av 

V 


Earth 


( 1 +Ae) - (Bl * e-e 1 2 ’.?) 

d+Ap') - itz't) 


ll i (1-Ap* ) 1 1 

Y 2 " 2 (1-Ae) 2 


(B4) 


Y3 = Y 1 r since we have assumed that the earth moves in a circle 
with a constant velocity. Expanding the denominator to fourth- 


order yields 


Av 

V 


Earth 


i 


[ (l-Ap'+Ae+Ap' ^-Ap' ^+Ap' ** -Ae Ap ' +AeAp ' ^-AeAx^' 

- (Bi*e) + (ei2*^) + (Bi’e)Ap' - (ei 2*^) Ap' 

“ (Bi*e) Ap'^ + (eiz*^) Ap'^ - (e*B) + (e*B) Ap' 

+ (^ 2 *?) - 2 (^ 2 *?) Ap' + (B 2 *e) Ae + 3 (B 2 *e) Ap'^ 

- 2 (e 2 *e) AeAp' - (B 2 *e) (Bi*e) + (Ba'S) (£ 12 *^) 

, Y3 

+ (32*e)^] x-j [2 - 1 + Ap ' - Ae + Ae Ap ' 

- Ae^ - Ae^ - Ae"* + Ap'Ae^ + Ap'Ae^] - Y • (^5) 
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The explicit functional form of Y3/Y2 can be taken from Eq. (A 4 ) 
of Appendix A; inserting this into Eq. (B 5 ) , multiplying out, and 
cancelling and combining terms, finally yields. 


V 


Earth 


= + I (31^-82") + Ae(Ap'-Ae) 

+ j [ (ei2*^) + (^2-?i) •£ ] 

+ 


- (Ap'-Ae) f ( 


m 


m 


-) + i ( 3 i^- 32 ^) + Ae (Ap*-Ae)] 


ri r^' 2 

~ 7 [ (£i 2*^) Ae + ("^2’e) Ap' - Ae + (c‘t) ] 

1 


+ ra [ y (-£J + 


3 3 1 1 32^ 


32^^ 


8 


1 2 

_ D ** 


2 r 


2 ^2 


[ 33T - 32“ - 231^32^] 


+ m^ [ (3 -y) (■ 


rz 


r 1 


i-) ■. 

2 ^ V- 2 




1 

2 


3i^ 


)] 


r r 
1 2 


) ] 


+ (Ap'-Ae) { Ap' [ (~ - ^).+ j ( 3 i^- 32 ^ 

1 ^2 ^ 


+ Ae (Ap ' -Ae ) + Ae ^ } 

+ [ (1.2-J) + (?2-la)-?] [(^ - 7 -) + I 

r 1 ra ^ 

+ ^ { - (^i *e) [ Ae (2 Ap'-Ae)] + {^2*£)I 2 Ae (Ap'-Ae) 

+ Ap'^] + (ei2*^) [ Ae (2 Ap'-Ae)] 

+ (£•?) Ae + (^2*e) (?i2*^) + (^2*c) ($2“?i)*e ) 
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APPENDIX C 


The general relation for the doppler-cancelled signals received 

With 5. H fo + t ' and e'l 2 = - e 0 1 + e ' , we have 

Ap" E "^0 • £oi 

= (^ 2 “^') •• ( — £12 + £*) 

= - B •£ + ^ *e' + e 4' - e '4’ 

2 12 2 12 

= - Ap' + ^2*e' + “ s'*^' . (C2) 

Ae " = 1 1 • e 0 1 ' 

= -^i-si2+.^i*e' 

= - Ae + . (C3) 

Substituting Eqs. (C2) and (C3) into (Cl), and expanding to 
fourth-order, yields 

= (1-Ap') (1+Ae+Ae ^+Ae ^+Ae ** ) x 

Y2 , Y2 

[ y7 ~ 2 y7 * 

- AeAp ' +AeAp ' ^ -Ae Ap ' ^ 


Av 

V Probe 


at the probe-based station is 


Probe 


^ d-Ap' ) ^ 

(1-Ae) ^ Yi 
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- - 2 Ap' (^2*e') ~ 2 Ap' (eiz*^') 

+ Ae (^ 2 *e') + Ae (ei 2 *^') 

+ Ap' (^1 •£ ' ) 

+ 2 Ap' (t’-t*) + + 2 (^2,-e') (ei2*^') 

+ (ei2*^')^ + 3 Ap'^ (B^2*e') 

+ 3 Ap'^ (ei 2 *^') - Ae (e'»^') 

- 2 AeAp' (32*e') - 2 AeAp' (£ 12 *?') 

” (^i*e') (^2*e') " (l^i*e') (ei2*^*) 

- Ap ' 2 (^ 1 •£'))]- ~ . (C4) 
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The explicit form for Y2/Y1 ■» can be immediately picked up from 

Eq, (A 4 ) of Appendix A. To find Y2/Y0 explicitly, we use 
^ 

r 2 = ro + p : 

2 2 2 
ri = r2 * r2 = ro + 2 ro * P + P 


= ro^ ( 1 + 


2 ro»p , p' 


ro 


ro 


Therefore, similarly to Eq. (A 4 ) , we have 

— = 1 + (2 E_) + i (B 2 ^-Bo^) 

Yo r2 ro 2 

32 62 o 

+ m [ V (_^ - -^) + 3 82 


1 3 0 ^ 1 , 3 0 3 2 


r 2 


ro 


2 r: 


2 ri 


- ^ ( 


2 'r; 


+ 


ro 


+ ^ [ 362** - 3o“ - 2g2^3oM 


+ m^ [ (3 -y) (- 


ro r2 


1 .) L_i 

z' 2^2 2^2 r2ro'^ 


r 2 ro 

2 


:: 1 - 5L. [ «p _ p^ ^ 3 (£o_!_P) ] 

ro 2 2 „2 

r 0 0 0 


+ (3o-^') + I (6' ") . 


(C 6 ) 


The fourth-order terms in Eq. (C 5 ) do not give any additional 
corrections in the perturbation expansion. 

Inserting Eq. (C 6 ) and the_ explicit form for Y2/Y1 into 
Eq. (C 4 ) finally yields 


Av 


Probe 


,m_ _ m . 
'ri r2^ 



( 3 i^- 32 ^)] 


Ap' (Ap' -Ae) 


- j [ (ei2-?') + (^2-?i) •£’ ] 

+ (Ap'-Ae) .[(^ - ^) + i ( 3 i^- 32 ^) + Ap'(Ap'-Ae)] 
r 1 r 2 ^ 


•)] 


(C 5 ) 
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+ I ■{ [3 Ap'-2 Ae] [(ei 2 *^') + (32*e')] 

+ (Ae-2 Ap') (3i*e') - (e'*6') 

1 ill _ 1 fill 

2 ri 2 

+ I [ 3 32 " - 31** - 231 ^ 32 " ] 

-i-1 

r ira 

+ Ap' (Ap'-Ae) - Ap'^ } 

~ ^ { (ei 2 *^') [ 5 Ap' (Ap'-Ae) + 2 Ae^ 

+ ( 2 ^2“^i)*e' + ei2*^'] + (e'*^') Ap' 

+ (^ 2 *e') [ 5 Ap' (Ap'-Ae) + 2 Ae^ + (^ 2 -^ 1 ) *e'] 

~ ( 3 1 • e ' ) [ Ae (Ae - Ap ' ) + 2 Ap ' ^ ] 

->■ 

2 (Ap'-Ae) (-?--^)] 
ro^ 


+ ^ [(-P!_) - I (£^P)^ 

io 2ro^ ^ ^ 

- 2 (Ap'-Ae) (io*^') } 


* 2r/n\/l 1\|3 1 IX 

+ [ (3-y) ( ) + -o -K 

+ (Ap'-Ae) { Ae [(^ - ^) + y ( 61 ^- 62 ^ 

r 1 r 2 2 


+ ^ (^^^) - (^o-"^' ) > 

""o ro^ 


[ Y ( 


3 2 

^ 2 


3 2 
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CAPTIONS 


Figure 2a. 


Figure 2b. 


Figure 3a. 

Figure 3b. 
Figure 4a. 

Figure 4b. 


Distance from the sun (in AU) , versus period, for an 
earth-launched heliocentric probe, that has a per- 
helion of its motion at 0.29 AU. (Period = 0.52 

I 

years) 

Curve a is the total redshift for an earth-launched 

heliocentric probe that has a perhelion of its motion 

at 0.29 AU, plotted as a function of the probe's 

period. (Actually - Av /v is plotted.) The 

infinite signal-propagation-velocity case is shown 

here. Curve b is the third-order ( in 3 ) terms that 

contribute to the total redshift. Curve c is the 

(negative of the) purely gravitational fourth-order 

2 

contribution (i.e., the m term). Some of the points 
at the beginning and end of the orbit have been 
omitted for simplicity; they either drop below the 
scale of the figure, or have opposite sign to the 
general curve. 

Same as Figure 2a, with perhelion =0.5 AU. (Period 
= 0.65 years) 

Same as Figure 2b, with perhelion =0.5 AU. 

Same as Figure 2a, with perhelion =0.1 AU. (Period 

= 0.41 years) 

Same as Figure 2b, with perhelion =0.1 AU. 
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